Active microrheology in the continuum limit: can the macrorheology be recovered? 
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Active microrheology differs from its passive counterpart in that the probe is actively forced 
through the material, rather than allowed to diffuse. Unlike in passive microrheology, active forcing 
allows the material to be driven out of equilibrium, and its nonlinear response to be probed. However, 
this also renders inoperable the fluctuation-dissipation theorem used to justify passive microrheology. 
Here we explore a question at the heart of active microrheology: are its results consistent with 
macrorheology? We study a simple model material - a generalized Newtonian fluid, with a small 
but arbitrary shear-rate-dependent component - and derive a general expression for dissipation 
due to probe motion, which remarkably does not require the non-Newtonian flow to be solved. 
We demonstrate that the straightforward application of active microrheology gives results that are 
inconsistent with macrorheology, even when the probe is large enough for material to behave as a 
continuum, unless the forcing is gentle enough to probe only the linear response. Regardless, each 
technique encodes information about the material; if suitably interpreted, the (macro-) constitutive 
relation can indeed be recovered from the microrheological data. We emphasize that more, rather 
than less, information would be obtained if the two methods disagree. 



For decades, conventional (macroscopic) rheology has 
been employed to characterize and understand complex 
fluids and materials Q. Some materials, however, would 
be impossible to study using macro-rheometry, due to 
the difficulty or expense in obtaining sufficient quantities 
for analysis. For such materials, techniques in microrhe- 
ology are being developed 0- 0- Q Most microrheol- 
ogy is passive, and relates the thermal fluctuations of 
translating 0, 0] or rotating 0, [^j colloidal probes 
to linear-response, near-equilibrium viscous and elastic 
moduli through the fluctuation-dissipation theorem. 

More recently, techniques in active and nonlinear mi- 
crorheology have begun to be explored. A colloidal probe 
is optically or magnetically pulled 0, 0, 0, Q, 0] or 
rotated [lfj through the material, and a 'microviscosity' 



F* 



n 



or r\ 



i 3 n* 



(1) 



is inferred, analogous to the generalized Stokes-Einstein 
relation assumed in passive microrheology. (Starred vari- 
ables will here refer to microrheological quanitites.) 

Active microrheology differs from its passive counter- 
part in it allows the material microstructure to be driven 
significantly out of equilibrium, providing a window into 
the material's nonlinear response. Such information is 
not available to passive microrheology, which by nature 
probes the near-equilibrium properties of the material. 
For example, active microrheology is expected to reveal 
shear/force-thinning and thickening in colloidal suspen- 
sions, as is seen macrorheologically (lilllij. 

A central question for microrheology is whether or not 
its results are consistent with those obtained using con- 
ventional macro-rheology. Passive microrheology seems 
to be on a firmer theoretical footing, as the fluctuation- 
dissipation theorem can be invoked 0]. Unfortunately, 
no corresponding theorem exists to ensure agreement be- 



tween macrorheology and active microrheology. While 
such agreement would certainly be reassuring, it is not, 
strictly speaking, necessary. In fact, we argue the oppo- 
site to be true. Any disagreement between the two tech- 
niques would yield more, rather than less, information 
about the material: each probes the material in a differ- 
ent way, and thus encodes different information. While 
further work would then be required to properly inter- 
pret microrheological data, any discrepancies encode ad- 
ditional physics which could, in principle, be extracted. 

Differences between active microrheology and its 
macro-rheological analog have already been noted. Stud- 
ies have thus far employed colloidal suspensions, as they 
are perhaps the cleanest and simplest systems for theo- 
retical and experimental study. For example, we demon- 
strated semi-quantitative agreement between micro- and 
macro-rheological measurements of the Brownian (mi- 
crostructural) component of the viscosity of a suspension 
|l3l Il4j . This agreement carried, however, an impor- 
tant caveat: the data must first be scaled (with colloidal 
volume fraction <f> for micro vs. with <fi 2 for macro) for 
a meaningful comparison to be made. Such agreement 
could only be obtained, however, because the proper scal- 
ing and physics for colloidal suspensions was known; for 
more complex (or unknown) materials, one would not 
know how to proceed a priori. 

Another significant difference involves non-continuum 
effects, which arise when the size of the probe is compara- 
ble to microstructural features of the material. This has 
been noted in falling-ball rheometry [T^ | and in active 
microrheology [lflj . and is not terribly suprising, since 
the material is not being probed in the continuum limit. 
Two-point microrheology was developed in order to ad- 
dress this issue in the passive microrheological context 
0, 01; a similar strategy may be useful as well for ac- 
tive microrheology More counter-intuitive, perhaps, is 
that the inferred viscosity depends on whether the probe 



2 



is pulled at constant force or constant velocity [l7l. fo r 
which a simple physical argument has been offered (13| . 

Here we approach the macro/micro question from the 
opposite standpoint: even when the probe is so large that 
the material behaves as a continuum, does microrhcol- 
ogy reproduce macrorheology? One might assume that 
it would, since one is probing the continuum properties of 
the material. On the other hand, the flow and stress fields 
established in active microrheology arc not viscomctric, 
but decay with distance from the probe. By contrast, 
(macro-) rheologists have long emphasized and utilized 
viscometric flows 0|- I n light of this distinction, any 
macro/micro agreement should appear surprising. 

To address this question in a simple, clear, and general 
fashion, we consider a probe that is driven through a 
generalized Newtonian fluid, whose viscosity 



i] = 7/0 + £771(7), 



(2) 



is nearly Newtonian, but which contains a small, shear- 
dependent component 771, whose functional form is arbi- 
trary (although by construction 771(0) = 0). 

In what follows, we derive a general relation for the 
dissipation due to a probe moving steadily in such a ma- 
terial. This allows the microviscosity 771 to be obtained 
without solving for the non-Newtonian component of the 
flow field. We demonstrate that the micro- and macro- 
viscosities are in fact different, so that a naive applica- 
tion of active microrheology via Q would be inappropri- 
ate unless the forcing was extremely weak. Instead, we 
show that 771 is related to 771 through a weighted average 
that depends on the shape and properties of the probe. 
This relation can be inverted numerically in general, and 
exactly in at least one case. 

We pose an expansion for the fluid velocity u — Uq + 
eu\ + where the leading-order solution simply solves 
Stokes' equations 

V ■ cr = 77oV 2 it - Vp = 0, V ■ m = (3) 
Uo(T) = U, u {r -> 00) -> (4) 

subject to no-slip boundary conditions. The 0(e) non- 
Newtonian correction then obeys 



V • <Ti = 77 V 2 iti — Vpi = — fi, V ■ Mi 
ui(r) = 0, u±(r — > 00) - 



= 
0. 



(5) 
(6) 



Here cr„ = —p n 8 + ?7o e ri is the stress tensor associated 
with the n th -order velocity field, e„ = Vm„ + (Vm„) t is 
the n th -order rate of strain tensor, and fi represents the 
body force fi = V • (771 e ). 

One can show that the dissipation due to a body mov- 
ing through a generalized Newtonian fluid is given by 



F*-U*= l -J 7?( 7 )e --edV^^J 7 7 ( 7 ) 7 W. 



Using J5J, the dissipation is given to 0(e) by 

F *' U * = \ j ( r ?o + e? ?i) e o : e Q dV + e J mei : e Q dV. (8) 

We now show the second integral to vanish. Writing 
Vo e i = °"i + P18, we note that the isotropic (pi) term 
vanishes upon contraction with (traceless) e , giving 
J 77 ei : e dV = J &\ : e dV. Expanding e and us- 
ing the divergence theorem gives 

J 77 ei : e dV = 2 J n ■ CTi ■ u dA + 2 J u ■ fidV, (9) 

where the surface integral is taken over the body and 
at infinity. Finally, a generalized reciprocal theorem for 
Stokes flow that accounts for body forces V ■ cr = — f , 



h-cri-UodA+ I Uo-fidV = / n-er -uidA+ \ Ui-i^dV, 

(10) 

shows © to vanish identically, since fo = everywhere 
and Mi = on all surfaces. We thus arrive at 

F* ■ U* = X - J (770 + e77i)e : e dV + 0(e 2 ) . (11) 

An analogous relation holds for rotating probes. Signifi- 
cantly, only the Stokes flow around the probe is required 
to obtain the dissipation to 0(e). 

A straightforward application of would yield 



e?7i = 



if - Vo = c f Vi(io) jjdV 
Vo Vo J HdV 



(12) 



Even in this simple example, the micro- and macro- 
viscosities are not equal unless 771 is nearly constant over 
the integration volume - i.e. 7 is very small (linear 
response, as with Brownian probes) or asymptotically 
large. Rather, 77^ represents a weighted average of rji, 
reflecting the non-viscometric flow around the probe. 

From now on, we scale lengths by the probe size a, 771 
by V01 7o by the maximum (Stokes) shear rate 7max , and 
denote dimensionless variables with tildes. We consider 7 
in 771 to be parametrized by 7c , and denote 5 = 7max /7c- 
With this notation, Eq. 1 121 becomes 



Ivi(S~io)^dV 



(13) 



We interpret the microviscosity in terms of 8, the (scaled) 
maximum shear rate from the Stokes flow solution. 

An obvious question is whether one can recover 771 
uniquely from 771. If not, then different macroviscosi- 
ties could yield the same micro- viscosity 771 . This would 
imply in turn that their difference 770 would obey 



(7) 



fl D {8~^)~ildV = {). 



(14) 
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n 


Translating 


Translating 


Rotating 






Sphere 


Bubble 


Sphere 


1 


7max 


3U/V2a 


V6U/a 




0.8 


Ax 


0.210 


0.250 


0.442 


0.6 


A 2 


0.092 


0.120 


0.267 


A 3 


0.056 


0.071 


0.184 


0.4 


Ai 


0.037 


0.048 


0.137 


0.2 


A 5 


0.027 


0.034 


0.107 





TABLE I: Prefactors (from eq. [T^J for the power-law viscos- 
ity fjl obtained active microrheology, given a macro-viscosity 
fj 1 ~ j n . Three different methods of forcing are presented: a 
translating sphere and bubble, and a rotating sphere. 



Because fju is a function of S and 7 2 is positive definite, 
all coefficients of a power series in <5 of f)o must be zero 
for this relation to hold. This would imply fju = 0, and 
we conclude that l|13|) is indeed invertiblc. 

Having arrived at the general relation, we explore il- 
lustrative examples of consitutive relations fji and probe 
forcing. The simplest example, a fluid with a slight 
power-law viscosity, is recovered microrheologically as 



rji = d 7 



A„ <f\ where A„ 



Ii 2 o +n dV 



(15) 



The functional form for f)i(5) is identical to 7/1(7), but 
with a probe-dependent prefactor A n . For example, a 
translating sphere (TS) has shear rate 



1 + (2 - 6f 2 + 3f 4 ) cos 2 9 



7max 



giving coefficients A n shown in Table Coefficients for 
translating bubble (TB), 



7 = f 2 cos6>, 7max = v / 6?7/a 
and rotating sphere (RS) probes 

7 = r~ 3 sin6>, 7 max = 3V2fJ. 



(17) 



(18) 



are also shown in Table |U 

Even for the simplest model fluid, A n ^ 1, and so the 
microviscosity differs from the macroviscosity. Further- 
more, general viscosities fji can be expanded as power 
series in the low-<5 limit. 



E 

n=l 



B n 5 n r, (19) 

whereas an expansion of the microviscosity would give 

00 




FIG. 1: Microviscosities obtained from a fluid with abruptly 
shear-thickening viscosity 771 using a rotating sphere (RS), 
translating sphere (TS) and translating bubble (TB). Al- 
though the shear-thickening is seen to begin at the correct 
shear rate 8, the microviscosities approach the macro- viscosity 
only slowly. 



Although f)1 ^ fji, the macro- viscosity coefficients B n 
could be determined microrheologically, since A n can 
be considered known. Again, the straightforward ap- 
plication of active microrheology is inconsistent with 
macrorhcology, except in the linear-response limit (7 — > 

0) . Nonetheless, the connection between the two exists 
and can be extracted with proper interpretation. 

We now turn to the opposite limit, when the probe is 
pulled strongly enough to sample interesting features in 

1) 1 . The simplest illustrative function for 7)1 is a fluid that 
thickens abruptly, via a step function, 



m = e^-i/s). 



(21) 



For 6 < 1, nowhere does the material thicken, and one 
recovers fji (S < 1) = 0. When 5 > 1, some volume of 
the material thickens, as seen in Fig. ^ Eq. (|13|) can be 
evaluated exactly for TB and RS forcing, 



= 1 - 







>h 



5^/2 

[A5 2 - 1)VF 



1 

5^3' 



(22) 



1-36 2 



1 [Vs^T/d] 



45 3 



23) 



Although all three probes do show the onset of shear 
thickening at the correct S, they approach the 'true' value 
only slowly: like J™ 1 / 2 for translating bodies and S^ 1 
for rotating. The reason for this slow decay is purely 
geometric: the fluid thickens until 7 ~ 1/5, which occurs 
at r ~ 8 1 / 2 (r ~ (5 1 / 3 ) for translating (rotating) probes. 
The integrand scales like <5 -2 , and the volume element 
scales like S 3 ^ 2 (S) for translations (rotations), leaving a 
<5 -1 / 2 (<5 _1 ) decay. Unfortunately, such slow decays are 
ubiquitous in active microrheology, and are remnants of 
interesting features at intermediate values of 5. 

Rotating probes are marginally better at reproducing 
the macro- viscosity because the faster spatial decay of 7 
effectively weights the average closer to the probe. Tak- 
ing this strategy to its logical extreme (a shear rate that 



4 



decays infinitely quickly) would yield a microviscosity 
that averages strictly over the probe surface, and does 
not suffer from these long-range geometric effects. 

In fact, we demonstrate that this idea can actually 
be realized. Two of the probe motions described above 
(TB and RS) are special, in that the flow for different- 
sized probes are identical if their forces/torques are the 
same (meaning aU a = bUb and a 3 Q a = b 3 fl(,). Taking 
the difference between dissipation integrals in at the 
appropriate 5 a and 5b yields a volume integral over the 
difference in probe volumes, since the integrands arc iden- 
tical outside of the larger probe. For example, using this 
technique for the TB mode gives 



where 8„ 



(b 2 /a 2 )Sb to enforce F* 



F* 



(24) 
F*. The 



limit a — > b yields an expression for the surface-averaged 
macroviscosity, 



We close with a brief reflection on the agreement 
between (active) microrheology and macrorhcology, by 
posing a broader question about constitutive relations 
(whose validity we have assumed without comment). 
Such relations generally assume viscometric flow and 
large-amplitude, steady deformation 0|. It is possible 
that probes in active microrheology do not incur strains 
that are significant or steady enough for such constitutive 
relations to hold, despite the magnitude of 7. While the 
relations derived above may be confidence-inspiring, we 
feel it would be more interesting if they were experimen- 
tally shown not to hold. In that case, active microrheol- 
ogy would in fact be providing information that is inac- 
cessible to macrorhcological techniques. Since most flows 
in the real world are not viscometric, it would be inter- 
esting and important to investigate what accuracy (or 
relevance) constitutive relations hold in more general sit- 
uations. Active microrheology seems ideally suited for 
such investigations. 



25 



df,i 



05 



f,t TB = 



4tt ./ f= i 



where 6b = 5. Since 5j — 5 cos 9 = 
perform the integral over <fi to obtain 



(25) 



1, we 



~*TB 



25 



OS 



vl TB 



s 3 



fji{u)u 2 du, (26) 



which can be inverted exactly to reveal a remarkable re- 
lation between the micro- and macro viscosities, 



)_d_ 

S 2 86 



(27) 



One can verify l|27|) for the step-function and power- 
law viscosities considered above. The analogous surface- 
averaged viscosity relation for RS forcing, 



d(iH RS 5) _ 3 
85 ~ 2 



771 (5 sin0) sin 3 6d9, 



(28) 



docs not appear to have such a simple inversion. 

Several important points emerge from this simple anal- 
ysis. Most significantly, the straightforward application 
of active microrheology will not yield results consistent 
with macrorheology, unless the applied force/torque are 
small enough that the material microstructure is not sig- 
nificantly deformed. In the nonlinear regime, the micro- 
and macro viscosities for weakly shear-dependent fluids 
are related through weighted averages that can be in- 
verted uniquely. This result holds only for weakly non- 
Newtonian fluids, whose velocity field is approximately 
Stokes flow, and whose non-Newtonian component is a 
small perturbation. Whether similar results exist for gen- 
eral materials remains an open question. 
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